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Abstract 

In this paper we present a statistical study of edge fluctuations taken with the Gas Puffing 
Imaging (GPI) diagnostics. We carry out a comparison of GPI signal from an extensive database 
including four devices (two Tokamaks and two Reversed Field Pinches). The data are analyzed in 
terms of their statistical moments Skewness and Kurtosis, as done in [B. Labit, et al, Phys. Rev. 
Lett. 98, 255002 (2007)]. The data align along parabolic curves, although different from machine 
to machine, with some spread around the best-fitting curve. A discussion about the meaning of 
the parabolic trend as well as the departure of real data from it is provided. A phenomenological 
model is finally provided, attempting to accomodate experimental evidence. 

PACS numbers: 52.35.Ra, 52.25.Xz, 52.35.Kt, 52.55.Hc 
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I. INTRODUCTION 



Edge plasma turbulence represents a main obstacle towards a good confinement of mag- 
netic fusion plasmas, and has been therefore studied for several years. Its understanding 
has considerably progressed, and nowadays the results produced by first-principle drift- 
interchange turbulence numerical models share several statistical features with experimental 

n 

ones PJ . However, we have not yet grasped a complete comprehension of this problem. 
One feature of edge turbulence phenomenology that has arisen interest of researchers is 
its (supposed) universality: statistical properties of edge turbulence appear quite similar 
regardless of the device they are measured in. Indeed, if the driving of turbulence is every- 
where provided by the same instabilities, it appears quite natural that some affinities must 
appear. Attempts of applying principles of universality to edge plasmas may be traced back 
since the works on Self-Organized-Criticality ideas 2]- It was however soon realized that, if 
one goes beyond the qualitative similarity of the signal, at the quantitative level relatively 
large differences apparently exist between different machines fsee, e.g. the lack of consensus 
about the analytical form of the PDF of fluctuations) Differences are 

unavoidable, since the different devices span widely differing ranges of operating conditions 
(geometry, magnetic field, temperature, ...). However, by focussing on those features that 
are common to all experiments, one may hope of disentangling the physical mechanisms 
driving the turbulence from the contingent details, related to specific geometrical and/or 
operational features of each machine. 

One purpose of the present paper is to present an extensive comparison of edge turbulence 
data taken from machines in different configurations: two Tokamaks (NSTX and Alcator 



C-Mod 



10|) and two Reversed Field Pinches (TPE-RX and RFX-mod [ij). The data 



were taken by our group using the same technique: Gas Puffing Imaging (GPI). This allows 
a straightforward intermachine comparison. 

The most informative way of expressing experimental time series is in terms of their probabil- 
ity density function (PDF). PDFs from all devices do indeed share some broad features: by 
considering density fluctuation PDFs, they are strongly skewed curves with approximately 
exponential tails towards high values, definitely departing from the Gaussian paradigm of 
independently distributed small fluctuations. When dealing with large databases, however, 
using the whole PDF becomes cumbersome, and more compact ways of expressing informa- 
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tion must be found. Quite recently, Labit et al [6| (hereafter referred to as Labit for short), in 



an analysis of TORPEX 



13| data, proposed of expressing data in terms of their normalized 



third (Skewness S) and fourth statistical moment (Flatness F). Besides the convenience of 
replacing a whole time series made of thousands of figures with a single couple of numbers, 
it was there argued that S and F do actually contain a large part of the physics embedded 
in the raw data in terms of their mutual relation: indeed, empirical data from TORPEX as 
well as from many other physical systems [l^ , together with general analytical results, point 
to the fact that S and F are roughly related by a quadratic relation F = A-S'^ + B- S + C, 
and the details about the physics driving the system are packed into the coefficients A, 5, C 
(although the linear term B is usually discarded, either because it is empirically found neg- 
ligible, or because theoretical considerations suggest that it should be null. This issue will 
be more extensively discussed later). 

In order to provide a physical rationale to the behaviour of a raw time series, one can be con- 
tent with just a phenomenological approach, i.e., picking up an analytical approximation to 
the empirical PDF, on the basis of the good agreement with data, but also backed up by some 
qualitative modelling of the physics running in background. This is the approach followed 
in the Labit paper as well as in several earlier studies [s], 0, [s, 3]- The coefficients A, B, C 
may be compared with those directly computed from the guessed PDF. Alternatively, a 
first-principle approach may also be attempted, where some more fundamental equations 
are postulated to model the system studied. Time series may be generated by solving these 
equations, and the corresponding statistical moments evaluated. Fairly recently, Krommes 
3. adapting to edge plasma turbu.enee an earlier model originally developed to study 
the fluctuations of sea surface temperature [1^ was able to reproduce a F-versus-S curve 
closely resembling TORPEX data 6|. It is particularly remarkable that, although the rele- 
vant equations are stochastic nonlinear differential equations, the computation of S and F 
in terms of the parameters of the model may be done in analytical closed form. 
In this work the data will be analyzed in terms of the moments S, F and results will be 
compared with Labit'. We anticipate here that some differences with those works will be 
found and discussed. The last part of the paper will attempt to organize these novel features 
within a phenomenological model previously advanced j^. 
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II. THE GPI DIAGNOSTICS 



Before entering the discussion of the results, let us provide some more details about the 
instrumental technique and data-taking procedure (more detailed reports are found in IgI). 
The GPI diagnostic is an optical, non-intrusive diagnostic for studying the plasma edge of 
fusion devices, and particularly the edge turbulence. It measures the visible light emission, 
usually emitted by neutral Hydrogen, Deuterium or Helium puffed into the plasma edge. 



In the two RF 
scribed in ref. 



experiments (RFX-mod and TPE-RX) the same GPI has been used, de- 



171]. It consists of a gas-puffing nozzle, and an optical system for collecting 



the light emission. The optical diagnostic arrangement views edge regions in the plane per- 
pendicular to the main component of the magnetic field, that is the radial-toroidal plane in 
the RFP experiments. It covers an observation region of about 80 mm along the toroidal 
direction and 40 mm in the radial one, by means of 32 lines of sight (LoS). The signals 
are sampled at 10 Msamples/s. GPI data are line-integrated measurements, although the 
effective length of integration is limited by the volume where emission of the puffed neutral 
gas takes place, which is very small (order a few cm). We note that this might a difference 
with respect to electrostatic probes, which sample smaller regions of space. This issue will 
be the subject of further investigations. The two-dimensional mapping of the signal is ob- 
tained through a tomographic inversion 18|. In this work, however, we will present just the 
signal from one chord. In order to be sure that conclusions do not depend upon the chord 
chosen, we made the following test: since we are interested in the third and fourth moment 
(skewness and flatness) of timeseries, we computed them for all the RFX-mod pulses and 
for the seven central chords, and estimated for each septet the mean value < 5* >, < F >, 
as well as the dispersion around this value. It turns out that 95% of points have a relative 
dispersion in S lesser than 32% and a dispersion in F lesser than 18%. 

In NSTX device the GPI consists of a double array of photomultiplier tubes that observes 
the radial-poloidal plane, i.e., the plane perpendicular to the magnetic field at the edge of 



a tokamak device (see references 



IS 



20[|). The signals are sampled at 0.5 Msamples/s. The 



GPI installed in Alcator C-Mod is quite similar, since it has 2 optical arrays that view the 
local gas puff in the plane perpendicular to the main magnetic field, and the light is collected 



by photodiodes 



211 ]. The sampling frequency is 1 Msample/s. 



The different arrangements imply that RFX-mod and TPE-RX measurements are line inte- 
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grals along the radial direction (with an effective width of order 2-3 cm), while NSTX and 
Alcator ones are along the toroidal direction. In principle, this could affect our results if 
measured quantities were strongly varying along the radial direction over distances smaller 
than about 1 cm. 

The whole database consists for RFX-mod of 2167 plasma discharges with different plasma 
currents and densities; for TPE-RX, 40 plasma discharges at low plasma current (150- 
350 kA) have been considered; for NSTX only 13 discharges with plasma current of about 
800 kA and with a spontaneous L-H transition have been used. In this way, both the L- and 
the H-mode could be studied in the same set of plasma shots. All the available GPI signals 
of NSTX have been considered, and so both the poloidal and radial array of chords have 
been analyzed. The data are collected in different radial positions, from about 60 mm inside 
the separatrix up to about 60 mm outside. Finally, we collected data from 10 Alcator dis- 
charges. Regarding typical plasma conditions, electron temperature in the collection volume 
is in the range 10-40 eV for all devices, while typical plasma density is of order 10^^ m~^ for 
RFX-mod, Alcator and NSTX, and up to an order of magnitude smaller for TPE-RX. All 
timeseries were splitted in strips 10 milliseconds long, and the couple {S, F) was computed 
for each strip. Since the available part of the discharges is usually much longer than 10 ms, 
more than one {S, F) couple was computed for each discharge. The time interval was chosen 
as a good compromise between competing demands: it had to be as small as possible in 
order to consider as almost stationary the background plasma conditions during the single 
time series, but long enough to prevent undesired fluctuations due just to the scarce number 
of points, as well as spurious correlations due to microturbulence characteristic times, as 
guessed from power spectrum. Depending upon the device, hence, each couple {S, F) was 
computed using up to 10^ points (for RFX-mod) down to 5000 (for NSTX). For comparison 
purposes, some trials were done, using RFX-mod data, by splitting time series into strips 
only 2 ms long. No relevant differerences were spotted on statistical terms: the average nor- 
malized difference < \SiQms — S2ms\/ Sioms > turns out to be less than 20% (the subscripts 
refer to moments computed using respectively the longer and the shorter time strips, and 
the brackets <> imply an average over all pulses). 
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III. INTERMACHINE COMPARISON 



The main result of this section is summarized in the figures below, where we plot the 
database of couples {S,F). Each subplot contains the database from one device. Overplot- 
ted, there is a subset of the RFX-mod database selected so as to keep matched a few oper- 
ation parameters in both devices. This way, we attempted to compare data from roughly 
similar experimental conditions. In RFPs three meaningful parameters are: (I) the ratio 
N/Ip between the (line-integrated) density and the plasma current, or equivalently n/na 



between the average density n and the Greenwald density no [23|]; (II) the reversal param- 
eter F = Btor{o)/ < Btor >, ratio between the toroidal magnetic field at the edge and 
its value averaged over the whole plasma. In RFPs it is negative, since the toroidal field 
reverses sign close to the edge, and is used to parametrize the equilibrium. (Ill) Finally, 
the absolute value of plasma current Ip is an important parameter, too. TPE-RX data 
(and the accompanying RFX-mod ones) are parametrized by: (I) low normalized density, 
0.1 < n/nc < 0.3; (II) deep reversal parameter, F < —0.08; (III) low plasma current, Ip 
lower than 350 kA. In tokamaks, the reversal parameter F does not make sense since it is 
trivially fixed by the aspect ratio. Hence, NSTX pulses are characterized by the normalized 
density: 0.25 < n/ric < 0.45 and the plasma current Ip = 800 kA. Alcator C-Mod pulses, 
finally, have 0.1 < n/nc < 0.5, realized in correspondence of two values of Ip = 0.4,0.8 
MA (corresponding to two values of magnetic field: 2.4 and 5.4 T). In this case we did not 
overplot RFX data. 

Several conclusions may be drawn from the plots ([T]): data do indeed distribute around 
quadratic curves, as shown by Labit. However, there are sensible differences between ex- 
periments. Loosely speaking, results may be parametrized according to the degree of order 
that is commonly associated to each magnetic configuration: tokamaks data-with Alcator 
in primis-align fairly well along a unique curve; NSTX in L-mode does the same, while in 
H-mode there is slightly more scattering at large 5*. Interestingly, the L-mode and H-mode 
curves are definitely different. This is an hint that the present analysis is sensitive to changes 
in the plasma transport. RFPs data are unequivocably scattered throughout a finite area. 
TPE-RX data appear more concentrated, but it could be an artefact due to the poorer 
statistics. At this stage, we cannot still conclude whether the spread of the data around 
the interpolating curve is a consequence of the non-homogeneous set of plasma conditions. 
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FIG. 1: (Color online) Top-left panel: yellow triangles, TPE data; dots, RFX-mod data selected 
so as to match TPE parameters (see main text). Solid curves are interpolating second-degree 
polynomials F = A-S'^^B-S^C^ and the values of the coefficients A are explicit displayed. 
Top-right panel, the same for NSTX. Here, both the data in L-mode (yellow) and H-mode (red) 
are displayed. RFX-mod data are selected to match NSTX parameters. Bottom panel, Alcator 
C-Mod data with interpolating polynomial. 



although it appears a plausible conclusion. 

In order to make more quantitative estimates, we have tabulated into tables I and II the 
numerical values of the fitting coefficients. In table I we reported the results for the most 
general fitting: F = ylxS'^ + -BxS' + C, while in Table II we did the same but postulating 
that the linear term were not present, just like in Labit: F = A' x S*^ -|- C . 

Obviously, the three-parameter fitting performs better in the least-squares sense than 
the two-parameters one. However, it is interesting to notice that even the three-parameters 
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experiment 


A 


B 


C 


RFX-mod (plot a) 


1.27 ±0.05 


0.91 ±0.11 


2.44 ± 0.06 


TPE-RX 


1.14 ±0.20 


2.03 ±0.66 


1.17 ±0.52 


RFX-mod (plot b) 


1.29 ±0.01 


0.82 ±0.03 


2.57 ±0.01 


NSTX (L-mode) 


1.74 ±0.31 


-0.98 ±0.77 


3.42 ± 0.46 


NSTX (H-mode) 


1.13 ±0.40 


2.75 ±0.93 


1.89 ±0.48 


Alcator 


1.48 ±0.04 


0.35 ±0.07 


2.88 ±0.03 



TABLE I: Fitting coefficients for the full parabolic fit F = Ay. S"^ + B y. S + C . The rows RFX-mod 
(plot a) and RFX-mod (plot b) refer to the two subsets of RFX-mod data compared respectively 
to TPE-RX and NSTX data. 



experiment 


A' 


C 


RFX-mod (plot a) 


1.65 ±0.02 


2.90 ±0.02 


TPE-RX 


1.75 ±0.04 


2.73 ±0.11 


RFX-mod (plot b) 


1.67 ±0.005 


2.92 ± 0.005 


NSTX (L-mode) 


1.35 ±0.04 


2.85 ±0.08 


NSTX (H-mode) 


2.28 ±0.10 


3.19 ±0.20 


Alcator 


1.62 ±0.02 


2.97 ±0.02 



TABLE IL Fitting coefficients just hke in table I, for the reduced parabolic fit F = A' x S*^ ± C' . 



fitting is not able to catch the whole physics embedded in the data: if data do really align 
along a unique parabolic curve, then the residuals would be of purely statistical origin, and 
therefore would be distributed according to a normal distribution. D'agostino-Pearson's is a 
statistical test designed to assess whether a sample of data derives from a normal distribution 
jsoj l. We carried on such a test on our data and found that it yielded a negative answer (to 
within the standard 0.05 significance level): the scatter of data around the quadratic fitting 
curve is not distributed as expected from a normal distribution. We note that some degree 



of scattering is present also in TORPEX data 



6j but was not given a particular emphasis in 



that analysis. We recall that we are not dealing here with raw data but with moments, i.e., 
averages of the data. This means that fluctuations of purely statistical origin should-to a 
large extent-be washed out by the averaging procedure: it is easy to show that the sampling 
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errors on S, F due to the finite length of the time series are orders of magnitude smaller than 
the observed scatter. Therefore, there must be another cause for the spreading of {S,F), 
not related to insufficient statistics. An obvious guess is relating S, F to microscopic plasma 
conditions that are fluctuating from one measurement to another. 

Also remarkable is that the coefficient of the quadratic term is always close to, but usually 
significantly different from, the value A = 1.5 recovered for TORPEX data. The value A = 
3/2 was first justified in Labit arguing that the functional form of edge density fiuctuations 
were closely matched by a Beta distribution 2J]: an analytical curve often encountered 
in modelling statistics problems. The parabolic relation between S and F, by itself, is not 
surprising; it is often reported in experiments from different environments (from atmospheric 
sciences to oceanography [l^). Indeed, speculations may be put forth, supporting the view 



that it must be a fairly general result. Very recently, Krommes [15| has developed a model for 
studying edge density fiuctuations. An exact result of him is that F = (3/2 — ao)5'^ + C(ao), 
with ao a free parameter within the model. 

Within Krommes' modelization, a change of F-versus-S trend may be taken into account 
through the coefficient ao, that potentially may vary between any couple of measurements. 
On the contrary, Labit analysis is rigid in prescribing A = 3/2. We further note that the 
interpolating curves feature non-negligible linear terms B ■ S. This is particularly evident 
for NSTX in H-mode, but is generically true for all the fits. Again, this result clashes with 
Labit' modelling using Beta distributions, but in this case also Krommes' model fails, since 
no linear term is there expected. 

NSTX fits feature striking differences between L-mode and H-mode, suggesting deep changes 
to the underlying physics. Table III points that, in H-mode, the {S, F) is acquiring a stronger 



experiment 


A' 


C 


A' 


B" 


C" 


A" 


NSTX (L-mode) 


1.35 ±0.04 


2.85 ± 0.08 


4.44 


3.31 ±0.13 


0.97 ±0.16 


6.72 


NSTX (H-mode) 


2.28 ±0.10 


3.19 ±0.20 


106 


5.28 ±0.25 


3.19 ±0.21 


105 



TABLE III: Fitting coefficients for the reduced parabolic fits F = A' x S'^ + C (second and third 
column) and F = B" x S + C" (fifth and sixth column). A' and A" are the sum of residuals 
squared (A = X]j((ij — c^ijit)^) for the two fits respectively. 



and stronger dependence from the linear term at the expenses of the quadratic one: the linear 



9 



fit perform as well as the quadratic one. 

A further interesting issue is that of the existence of time series with negative skewness. 
In our measurements they are very unlikely, and furthermore are rather marginal: —1<S 
(see Fig. 2). Physically, negative skewness in density data implies that fluctuations are 
favoured depleting particles from the region of the measurements. The possibility of their 
existence was documented by DIII-D 22 1, TORPEX 6], and Large Plasma Device 2^, that 
interpreted them as density holes. The exact reason of this discrepancy between our results 
and theirs is as yet unclear. We note that all of the other mentioned analysis have been 
carried on using edge Langmuir probes. The issue of the correspondence between quantities 
measured by probes and those measured by GPI is not still completely settled. Up to now, 
all investigations done by our group have not yet found anything invalidating the view that 



essentially the same physics is monitored by both diagnostics [26|. On the other hand, quite 
recently theoretical speculations have been put forth, suggesting that in some situations 
definite differences may be expected 

The above results suggest therefore that a fully satisfactory modelling of edge fluctuations, 
valid for all devices, has not yet been achieved. In the next section we reconsider an instance 
of such a modelization, advanced in an earlier work, and show to what extent it may fit the 
present results. 

IV. AN ATTEMPT OF PHENOMENOLOGICAL MODELLING OF EDGE FLUC- 
TUATIONS 

A first attempt of phenomenological modelling of GPI data in RFX-mod was given in . 
The need for a different modelling tool with respect to existing proposals was stimulated by 
the recognition in RFX-mod data of a novel feature, unnoticed in other experiments, i.e., 
the frequent presence of a double slope in the tail of the PDF 5|. We may now confirm 
that a similar feature is spotted also in other devices (Alcator C-Mod and TPE-RX). In [sl 
we justified this double slope by postulating the existence of two contributions to the total 
signal, hence that total PDF be a linear combination of two basis functions: 

P{X) = CiPiix) + C2P2ix) (1) 
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^'^""^ = r(iV-) """P (-^^•^^•^) ' r(^) = t^-' exp(-t) dt (2) 



In the original work [5| we used Gamma distributions as basis functions: 

\Ni /-oo 



Gamma distributions are fairly versatile, and may fit a large variety of empirical curves. 
They satisfy an important constraint: are defined only for x > 0, hence prevent density 
being negative. This was somewhat overlooked in earlier studies, where fluctuations around 
the average value x— < x > were considered, rather than the whole signal, missing therefore 
an important information. 

The parameter Nj plays the role of a number of degrees of freedom: it can be shown that Pj 
is the resulting distribution when our signal is made by the sum of Nj independent stochas- 
tic variables yl''\l = l,...,Nj, each of them being distributed according to a Boltzmann 
distribution: 

P{y) = -^exp(--^) (3) 



<y> \ <y> 

Since Ni ^ N2, the two contributions in ([T]) were thought as coming from two mechanisms 
living in abstract spaces of different dimensionality. When ^ 00, Gamma distributions 
approach Gaussians, hence it is natural to relate the high-A^ contribution to turbulence 
coming from uncorrelated small scale fluctuations, while the low-A^ one is modelling some 
correlation that diminishes the effective number of degrees of freedom. See 



for an earlier 



suggestion pointing to this kind of mechanism. Note that Eq. ([T]) refers to addition of 
probabilities, not of densities. In other terms, we are postulating that, at each time step, 
the signal is due entirely either to one mechanism or to the other, but not to both. This is 
clearly just a first approximation. 

It is obvious that the choice of the basis function is constrained only by the final accor- 
dance with experiment, hence we may inspect the effect of using different basis function. 
Accordingly, we thought of straightforwardly generalizing TORPEX analysis, using as basis 
functions Pi,P2 two Beta distributions: 



B{pj,qj){xf^ -x^Py^+1^-^ 



Pj{x) — P{x,pj,qj,Xj ,Xj J - (^)^ 

B{p,q) = [ tP-\l-ty-^dt (4) 
Jo 



[Note that we use the symbols (3j within the definition of Gamma PDFs ([2]). This turns out a 
bit unfortunate in the present context because of possible confusion with Beta distributions, 

11 



but is needed in order to be consistent with earlier literature]. 
From Eq. ([1]) we then evaluate the mean value and higher moments: 



< X > = J xP{x)dx 

M{i) = J{x- < X >yP{x)dx,i> 1 (5) 

in particular, 5" = M(3)/M(2)3/2, F = M{A)/M{2f. Although exphcit analytical formulas 
may be given for S and F when Pj are given by Eqns. ([2]) or (jlj), we won't write them 
down here, since they are fairly lengthy. It is however obvious that expressions ([5]) depend 
on a large number of free parameters: there are six of them in the case of two Gamma 
distributions, but two are fixed by requiring normalization of P: Ci + C2 = 1, and of average 
density: Ci//3i + C2//32 = 1- For Beta distributions, the situation is more intricate, since 
we start with a larger number of parameters: ten. We still get rid of one of them by the 
normalization condition on P. An interesting property of one single Beta distribution is that 
S*, -F, do not depend upon the boundaries x^''\x^^\ In the two-PDF case, this result does not 
hold any longer because of the way the dependences upon the parameters are hardwired into 
< X >. In order to simplify algebraic expressions, we considered < x > as a supplementary 
free parameter. Since it is formally independent of {x^^Xj'*^}, these latter quantities could 
therefore be discarded from final results. We are thus left with six free parameters. In 
principle, therefore. Beta distributions have more room for interpolation. 
Figure ([2]) summarizes the result of the analysis. Symbols are all the measurements from the 
three devices RFX-mod, TPE-RX and NSTX 29| , the shaded regions fill the parameter space 



spanned by the analytical curves ([T]) through a scan of the free parameters over their domains. 
It is interesting to notice that, although one Gamma PDF cannot have negative skewness, a 
linear combination of two Gamma PDFs can-at least marginally. Notwithstanding the fact 
that fitting with Beta PDFs ought be easier, since we have more free parameters available, 
actually these functions turn out to be less flexible: there is not a large difference between 
using just one or a sum of two Beta PDFs. 

On the whole, the blue region overlaps better the cloud of experimental data: most of them 
fall inside it. On the contrary, about one half of the points lies outside the orange region. 
Hence, we can claim that using a linear combination of two Gamma PDFs should, on the 
average, perform better on our data. However, the purpose of our work was not as much 
assessing which analytical PDF works better (in principle, there could be other better choices 
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FIG. 2: (Color online) Symbols, GPI data from different experiments: grey crosses, from RFX- 
mod; yellow squares, TPE-RX; red crossed, NSTX in H-mode; green crosses, NSTX in L-mode. 
The blue region fills the area spanned by the two-Gamma PDFs; orange region is the two-Beta 
PDFs. 

until now unexplored) as establishing the minimal degree of sophistication needed for a good 
modelling of edge fluctuations. 

V. SUMMARY AND CONCLUSIONS 

Summarizing, we have reached the following conclusions: (I) The physically relevant issue 
of the universality of edge turbulence between different devices has been addressed, with 
mixed results. From the one hand, we may confirm earlier studies: all data align along 
parabolic curves. However, we pointed out that the mere existence of a quadratic relation 
between S and F does not provide much insight. The relevant information is contained in 
the coefficients of the curve, that appear not univocal: they are very similar between two 
RFPs experiments (TPE-RX and RFX-mod), differ between a tokamak and a RFP (NSTX 
and RFX-mod), as well as between two modes of operation of the same tokamak (NSTX 
in L-mode and H-mode), and, finally between two tokamaks (NSTX and Alcator C-Mod). 
Furthermore, within one and the same machine, the spread of points in the (S*, F) plane 
warns that plasma parameters varying within the same shot strongly affect the statistics of 
the turbulence. (II) The analysis in terms of (5, F) introduced in Labit confirms to be a 
very effective and compact tool. (HI) As a consequence of (I), any attempt of modelling 
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phenomenologically the statistics of edge turbulence must be fairly flexible. The model 
advanced in the previous section appears promising under this respect. 
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